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SHORT METHOD OF ELLIPTIC FUNCTIONS. 



BY LEVI W. MEECH, A. M., HARTFORD, CONN. 

The system of Elliptic Functions is closely connected with analytical ge- 
ometry and trigonometry, and may be said to rank next in importance. The 
earlier advances of Fagnani, Euler, Landen and Lagrange were expanded in 
a systematic treatise by Legendre, whose well directed labor of forty years 
is also allied with the brilliant discoveries of Abel and Jacobi, and many 
others. Their researches have greatly facilitated the solution of numerous 
problems in Mechanics, Astronomy, Geometry, and the theory of Heat and 
Light. 

The present compendium has been arranged as free from complexity as 
possible. Especially in the seventh Section, the unexpected discovery of a 
new scale of amplitudes adapted equally to successive or to simultaneous 
quadrature, happily attains the object of the theorems of Jacobi and Ivory, 
without the previous long and intricate demonstration. And the common 
origin of all known scales of moduli is traced to a simple integral in the 
next Section. 

Section I. Change from Algebraic to Elliptic Functions. 

" II. Connection with the Spherical Triangle. 

" III. Increasing Scales of Moduli. 

" IV. Decreasing Scales of Moduli. 

" V. Integration by series. Table of Quadrants. 

" VI. Inverse Method of Integration. Theta Functions. 

" VII. Successive and Simultaneous Quadratures. 

" VIII. Quadrantal Quadrature. 

" IX. Elliptic Functions of the Second and Third Order. 
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Definitions. Elliptic Functions, or Transcendents, were classed by 
Legendre into three species. And the complete periphery of each is a clo- 
sed curve having four similar quadrants. The Function of the first species 
is denoted by F. The second, which is an arc of the common ellipse, is 
denoted by E; and the third by II. Thus, 

F=F(e,0)= P /n d \- mv 
J i/[l — e 2 sm 2 0] 

E = E{e, 6) = f 9 t/[1 — e 2 sin 2 0] dO. 

n = //(e, n, 0) = J o(1+nsin ^ )l/[1 _ e2sin ^ ] - 

To abbreviate, let A = i/(l— e 2 sin 2 0) ; 6 = -j/(l— e 2 ); i^e, J^) = A.\n. 

The variable arc is named the amplitude of the function. And e, being 
a fraction not greater than 1, is called the modulus, or, eccentricity ; while 6 
is termed the complementary modulus. Also n is termed the parameter. 
If e is denoted by sin <p , and 6 by cos p, then cp is sometimes called the an- 
gle of the modulus. 

Examples of F. From Legendre's Fonctions Elliptiques, Vol. 2, Table 9 : 
e. = 25°. = 45°. = 70°. = 90°. 



Sin 25° 

« 45 o 

" 70° 
" 90° 



25° 
45° 

70° 
90° 



0.43874 7923 
0.44328 23288 
0.44903 50293 
0.45087 533 



0.44339 4028 
0.42957 52476 
0.42242 64952 
0.42261 826 



0.79870 514 
0.82601 7876 
0.86652 9957 
0.88137 359 

Values of E. 

0.77247 109 
0.74818 6504 
0.71714 7672 
0.70710 6781 



1.26548 460 
1.36971 9477 
1.59590 6244 
1.73541 516 



1.18039 569 
1.09900 8292 
0.98297 5826 
0.93969 2621 



1.64899 522 
1.85407 4677 
2.50455 0079 



1.49811 493 
1.35064 3881 
1.11837 774 
1.00000 0000 



I. Change prom Algebraic to Elliptic Functions. This is gen- 
erally required by the occurrence of a radical denominator of the third or 
fourth degree. The first step is to resolve into two quadratic factors; thus 
the denominator B = 1 /[A+Bx+Cx 2 +Dx i + Ex*'] 
= i/[(a+6a;+ca; 2 )(a / +6'a!+c'a! 2 )]. 

The second step is to eliminate the odd powers of the variable. For this 
object, one method consists in changing to y 2 , after assuming 
a+bx-\-cx 2 = (a' -\-b'x-\-c'x i )y 2 . 

Another simpler assumption is, x — (p+qy)-+-(l-+y)- Substituting into 
the preceding factors, and equating each coefficient of y to zero, 
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a +lb {p+q)+cpq = 0, 

a'+%b'(p+q)+c'pq = 0. 
These give by elimination the sum p-\-q, and the product pq-, then^>, q 
are roots of the quadratic u 2 — (p-\-q)u-\-pq = 0. Should these roots come 
out imaginary, Legendre proves there are two other combinations of the bi- 
nomial factors (x — r) etc., of B, which would give real roots for p, q. The 
two factors of B thus take the required form f+gy 2 . 

1. When the factors of B, or (a? + 2ab cos s.x 2 + b 2 x*) H , are imaginary. 

Let x = ^(f)tan^, e = sin| S ; then g = ^y^^fy 

2. When B = \/[a\l -f 6 2 a; 2 )(l— eV)]. Let a; = (l-r-c)eos 6 ; 
= e ; which gives -^ — -t'- 



^.i^ *> 6 i? a&V[l— e 2 sin 2 0']' 

To make the result positive, let cot0 = ]/(l — e 2 )tan0'; then 

dx_ e dd> 

B ^VC 1 — e'sin 2 ^']' 

3. When B 2 = (1 +a?a?)(x t -b 2 ). Let x = * ; —1^ = e 2 ; ^ 

v ;v y cos0 l+a 2 6 2 jB 

~" t/[1— e 2 sin 2 <?]' 

4. When E 2 =(l +aV)(l +6W). Suppose a> 6 ; x=^- ; ^=^ = e 2 . 

ct a 

5. WhenJ? 2 =(l— aV)(l— 6V); a>6; let aa; == sin0; 6-r-a = e; then 

JE = _, — ^ _ „„, . The radical will be imaginary from a; = - to T , 

JS a 4/[l — e 2 sin 2 0] a b 

beyond which it becomes real to x = co. For these real values, B 2 = 

(aV— 1)(6V— 1). Let 6a; = -r— 3 ; - = e ; -^ = — ~ . -=. 

v A ' sin a B «V[1 — e 2 sin 2 0] 

6. When B 2 = (a 2 — a: 2 )(a; 2 — 6 2 ); a; lies between a and 5. If a > 6, 
a 2 — 6 2 62 efo 1 d0 



let e 2 p— , and a; 2 — — - -^.^ , ^ - a y j- 1 _ e 2 gin 2^- 

See Legendre's treatise, Vol. I, p. 4, .... p. 252. 

7. When B 2 = a+bx 2 +ox i +dx 6 ; let x 2 = z, etc. 

8. When B i = a+bxP+cx*; let i2 — z or xy, etc. 

9. When B? = a+6a;+ca: 2 +da; s ; make B = x.d y '+z, etc. 

10. When B 2 ==(a-)-&sm 2 0)(c+dsin 2 0); change first from sin0totan0. 

11. When B 2 = A+Bcosd+CsmO+DcofsW+Esmd cos(?+i^sin 2 0; 

make tan 10 — z, or sin0 = — - — 5 , cos0 = , , „ , d!0 = -— — 5 ; etc. 

* 1+Z 2 1+2* 1+2^ 
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OTHER INTEGRAL, FORMS. 

12. M 2 = {% — a)(x — h)(x — c) ; make one factor, as a; — a = y 2 . 

13, 14. f"*™ = ±,(F-E) ; f?*M = ^(E- Vf). 

1P . 1c AanW \ I ., » J\ rdd 1 „ e 2 sin0cos0 
15, 16. J ___= _^tanfl-^j ; J ^ = ^E- yj . 

17 Js^ = ^^_^_sin^ - 
., „ /"cos 2 # dd 1 / £7 et\ i sin cos # 

J A* 36 4 36 2 36M ^ ^ + JV* 

20 - Jot*- £("-'+«"-*)• 

21, 22. f~i = A tan d+F—E ■ (a tsm 2 0dd = A tan 6+F—2E. 

23. JA sinmO = — § A sin cos + ?^~- #+ ^ 2 .F. 

24. f#d» = JeM sin cos fl+ 2 + — E— *L F. 
J 6 6 

25. f . ^- g = 2J tan W+2F— 2.E, 
•/ Jcos 2 J0 

II. Connection with the Spherical Triangle. In spherical 
trigonometry it is proved that when one side and its opposite angle, as o, C, 
are regarded as constant (Chauvenet's Trigonometry, p. 239), we have 

da , db ,-. 

cos J. cos-B 

The constancy of c, O, makes the ratio e constant in the same triangle 5 
thus, sin-4 = esina; sin^B = esin&; sin (7= esinc; 

cos A=\/\l — e 2 sin 2 a], cobB=\/{1 — e 2 sin 2 6], cos C= — j/[l — ^sin 2 c]. 

If c denote the longest side, C will denote the greatest angle, and the 
equation below requires that cos C have a sign opposite to that of cos A and 
cos B. Now substituting the radical values in place of cos A and cos B 
above, and integrating, 

J y/\\— e 2 sin 2 o] + J i/[l— e 2 sin 2 6] = J C '^' 
For, had the three sides been regarded as variable, the differential equa- 
tion quoted above, would have a third term of the form C'de. As arbitra- 
ry constant here, it can contain neither a nor 6. When a=0 > b=-o; when 
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b = 0, a = c; in either case the last equation defines the form of the func- 
tion C ; thus, 

r da , C db_ C do 

J y\l—e*wa*ayJ i/[l— e 2 sin 2 6] J i/[l — e 2 sin 2 c]* 
Thus the function F of the first species has the remarkable property of 
being expressed by two similar terms, which are not separately integrable. 
The relation was first discovered by the celebrated Euler in 1761. 

In respect to E, or the second species, we differentiate the preceding equa- 
tion, clear of fractions and restore the cosines before described ; giving 

cos B cos Ccfa+cos A cos C.db +cos A cos B.dc = 0. 
Substituting for each product of cosines, its value from the well known for- 
mula of spherical trigonometry, 

cos A cos B = sin A sin B cos c — cos C, or 
= e 2 sin a sin b cos c — cos C; 
we find cos A.da-\-cos B.db-\~cos C.dc 

= e 2 [sin a sin b cos c.dc+sin a sin c cos b.db-\-sin b sin c cos a.dd] 
= e 2 d(sin a sin b sin c). 
Replacing the radical values of the cosines, and integrating, 

JVC 1 — e 2 sin 2 a]cfa-f JV[1— e 2 sin 2 6]c?6 = JV[1— e 2 sin 2 c]efc 

-j-e 2 sinosin6sinc. 
Or ■#(«)+ ■#(&) = i£(e)+ e 2 sin a sin 6 sin c. 

TAe equation of the Ellipse whose rectangular coordinates x, y, have their 
origin at the center, is 

VL 4- *! - 1 

J3 2 ^ J. 2 ~ 

Let denote an auxiliary angle ; and since sin 2 -f- cos 2 = 1, we can 
make y = B sin 0, a? = A cos ; B 2 — J. 2 (l— e 2 ). Then the arc 

l/[cfc 2 +<% 2 ] = Ai/[1— e 2 sin 2 0]d0. 
The latter member, when A is unity, is identical with dE; and the relation 
of the auxiliary arc 8 to the co-ordinates is manifest. Again if the prece- 
ding value of y is equated to that which is found from taking the focus near- 
est the apsis as origin, we find by comparison that 6 is the angle named in 
Astronomy the eccentric anomaly. And thus, if v denote the true anomaly, 
we have 

o!„ a l/(l— e*)sini> _ a e -f cos v 

sin a = -J^— ^ , cos a = - T — ■ ; 

1+ecosv 1+ecostt 

and the precding relation of three amplitudes or eccentric anomalies is ap- 
plicable to three distances on the orbit. 

The equation of the Hyperbola, whose rectangular co-ordinates have their 
origin at the center, is A 2 y 2 — B 2 x* = — A 2 B 2 , 
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If we make A-i-y'[A i -\-E i ']—e; y—B\/\l — e^tan 6; where 6 denotes 
an auxiliary angle ; x = j/[J. a -Hcos 2 — J. 2 e 2 tan a #] ; then the hyperbolic arc 

•BfaH-cy] = BVll-^ •^ di/ i f_ ehm = BVH-e^.dH. 

The Imaginary. Here let us change to 6', and after integration, return to 
; first assuming the symmetrical relation cos 6 cos 0' = 1. Then, sin d = 
!/(— I)tan0'; sintf' = y{— l)tan0; d0 = j/(— 1) ((20'h-cos 0') 5 

jjj. _ d d _ j/j—VcosW'dd' 

~ cos8 V[l — e 2 sin 2 0] j/[l — 6 2 sin 2 0']" 
Here 6 2 = 1 — e ! . By division and integration this becomes, 

H(e, 6) = ^j=±[E(b, d>)-e*F{b, 0')]. 

Eeplacing the symmetrical arcs 0, 0' by a, a', etc., adding the equations for 
a, b, and subtracting for c as the greatest side of the spherical triangle, the 
three terms of F become zero, as before shown, while those of E with accent- 
ed values change signs, and are equal to the product below ; 

H(e,a)+H{e, b) — H(e, o) = ^~ (5 2 sin a'sin&'sin c' W-tanatan&tanc. 

Taking e 2 = J, c = 70° and b = 45°, we find by the spheric formula be- 
low, a = 30°26'. Then H{6) = 3.2782, H(b) = 1.0617, H{a) = 0.6025; 
and each member becomes — 1.6140, which verifies the preceding formula. 
It therefore merits our confidence. The abstruse subject of integration 
through imaginary limits has been systematized by Cauchy and others. 

Among the earlier results, the following is readily proved by differentia- 
tion. 

TT— C dd -. ^tanfl — E+(l— e 2 )F 
J cos 2 0.J 1— e a 

To recapitulate, — when e is less than unity, let c denote the longest of 
the three sides a, b, c of the spherical triangle. Then will the three sides 
or amplitudes, and their elliptic integrals be connected by the following 
equations : 

cos c = cos a cos b — sin a sin &i/[l — e 2 sin 2 c], 
cos a = cos b cos c + sin 6 sin cj/[l — e 2 sin 2 a], 
cos 6 = cos a cos c -)- sin a sin cj/[l — e 2 sin 2 6]. 

F(a)+F(b) = F(c). 

E(a)+E(b) = i£(c)+e 2 sin a sin b sin c. 

H(a)+H(b) — H(e)— tan a tan 6 tan c. 

in \ , tt/i\ tti \ i n , -i r IV.sinasin&sinc 1 

11(a) + 11(b) = /7(c)+Tr=.tan - . r c . 

v ' ' w w N [_l-\-n{l — cos a cos 6 cose) J 
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Here N = j/[n(n+l)(m+e 2 )] ; and when n is negative, the last term is 
to be changed to a logarithmic form. Similar summations might evidently 
be obtained for the various integrals at the end of Section I. 

III. Increasing Scales op Moduli. The scale discovered by Lan- 

den, and extended to the limits by Lagrange, implicitly assumes that sin 

in one function F is proportional to the derivative of A in another func- 

i- r- t.\. i. ■ -e j. a sin 20' . a 2 sin 0'cos 0' 

tion i ; that is, it tan = — - ^ , sin r 



e + cos 20" 1 + e ' j/[l— e' 2 sin 2 0'] ' 

e' = h£L. Or sin(20'— 0) = esin0; e' 2 = 1— /t^\ 2 
1+e V ; U+«/ 

Or cos(20'— 0) = 4/p.— e 2 sin 2 0]; 2-5- (1 + e) *= e'-rVe. 
The agreement of the fourth and sixth of these equations is proved by 
squaring and adding ; and the agreement of the second equation with the 
first can be shown by squaring the first, then changing from sin to cot 0, 
and eliminating cot 0, which leaves an identical equation. 

Let us here introduce the auxiliary angle such that sin <p=e } sin <p'=e', . . . 
For computing the ascending scale ; 

p ,„Ve. „„_2iA' -,_2 1 /e". 

e -jz^' e -f+?' ~i+l"- ,m - 

or tanj^'-=i/e =-j/(sinf>); tan J{p"=j/(sin <p'); tan Jf"'=i/(sinf"); . . 
to be continued till the value of e is sufficiently near to the limit 1, which 
renders F an exact integral, as far as required. 
For computing the series of amplitudes, we have 

sin(20' — ) = e sin = sin <p sin ; 
sin(20" —0' ) = e'sin 0' = sin cp' sin 0> ; 
sin(20'"— 6") = e"sin 6"= sin y>" sin0". . . . 
Expanding sin (20' — d), and dividing the equation by sin 0, we find 
cotg= e+cos2g: ; dcrtg= -(l+eco sjng^ 



sin 20' * sin 2 20' 

XT 7 „ d0 — sin 2 0.c£cot0 — sin0.dcot0 

]N ow dB = — r = 



J J V[l + cot 2 0— e 2]' 

Substituting for sin its value from the assumed equation above, and for 
cot and d cot 0, the values just found, we have 

d F = dd == _? ^ 

^[l-AinVI l + eV[l— e /2 sin20'J 

When the scale e, e', e", .... has been extended so near to the limit 1 

that the last modulus may be taken as 1, let denote the corresponding 

amplitude, which is also a limit ; then 
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As the moduli increase, the amplitudes diminish to this limit ; thus, 

l?{e,0) = logtan(45 o +|0). 2 2 2 



Let if = 



l+e'l+e'l+e""* ' 
When = 90°, the first equation of amplitude reduces to sin (20' — 6}=l.e 
= sin <p, or 0' = 45°+Jf . 

_2 2 2_ - 1/ eVV" . . . 1 \ _ tan|p".tany'"...l 

1+el+e'l+e"" - ' AjV e / tanjp' 

JE[e, 0) = Jf.log tan (45° +|6>). 
The Napierian logarithm here indicated, is found by multiplying the 
common logarithm by 2.302585093. 

Preparatory to the next section, the difference of the two equations, sin 
= sin 0, and sin (20' — 0) = e sin 0, divided by their sum, gives tan (0 — 0') 
= [(1— e)H-(l+e)]tan0'= V[l— « /3 ]tan 0'. 

Note. The relation discovered by Gauss in 1818 gives the same moduli 
as before, but a different series of amplitudes; thus, e' = 2j/e-J-(l+e), 

gin >?/ - (!+*) sin fl . dO _ 1 <W_ 

l+esuM ' ^[1 — e s sin 2 0] 1 + e ' t/[1— e' 2 smW] 

This scale of Gauss can be derived from that of Lagrange by changing 
from and 0' to ¥ and ¥' after assuming cos cos ¥ = 1, cos 0' cos ¥'=1, 
as in Section II!. 

IV. Decreasing Scales op Moduli. The preceding scale of moduli 
can be prolonged in the opposite direction, tending to the limit 0; so that 
the entire scale will be 0, . . . e 000 , e 00 , e°, e, e', e", e'", ... 1. Therefore re- 
versing the accents in the above equation, we have for computing the de- 
creasing series of moduli, as well as the amplitudes: 

e — sin f, e° = sin y>°, . . .; e° — tan 2 J^>; e m = tan 2 J^°; . . . 

*. ao_2V&. 2 __ 6° 



^ = 1 /[l-6»»]=:±=^;6' 



1+6' ~ 1+6' 1+6 ~ Vb' 
tan(0° — ) — 6 tan =z cos f .tan , 
tan(0 00 — 0°) = 6° tan 0° — cos /.tan 0°, etc. 



*fc «) = ^. w fl") = ^°.^.^°°, n 



(To be continued.) 



